Stable Mean Field Solution of a Short-Range Interacting SO (3) Quantum Heisenberg 

Spin-Glass 
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We present a mean-field solution for a quantum, short-range interacting, disordered, S0(3) Heisen- 
berg spin model, in which the Gaussian distribution of couplings is centered in an AF coupling J > 0, 
and which, for weak disorder, can be treated as a perturbation of the pure AF Heisenberg system. 
The phase diagram contains, apart from a Neel phase at T = 0, spin-glass and paramagnetic phases 
whose thermodynamic stability is demonstrated by an analysis of the Hessian matrix of the free- 
energy. The magnetic susceptibilities exhibit the typical cusp of a spin-glass transition. 

PACS numbers: 75.50.Lk 



Spin glasses (SG) have attracted a lot of interest since a 
long time [l], Q • They are characterized by having some 
degree of frustration, which is caused by the competi- 
tion between different types of order. As a consequence 
some of their properties are shared with paramagnetic 
states and some other with ordered, ferromagnetic or 
Neel states. This is usually produced by a random distri- 
bution of coupling constants that allow for interactions 
of opposite signs. The characteristic time scale of this 
is typically much larger than the dynamical time scale 
intrinsic to the system, a fact that leads to the so-called 
"quenched" thermodynamical description [l], Q • 

In a pioneering work, Edwards and Anderson (EA) 
proposed a model for SG, introduced an order param- 
eter for the detection of a SG phase and employed the 
replica method to deal with the quenched average [3j. A 
simplified version of the EA model describing classical 
Ising spins with long-range interactions was solved for 
the first time by Sherrington and Kirkpatrick (SK) us- 
ing the mean-field approach 4}. The solution, however, 
proved to be unstable Q and this fact has been generally 
attributed to the so-called replica symmetry it possesses. 
Indeed, stable replica symmetry breaking solutions were 
subsequently found [|| and a lot of knowledge has been 
gathered about long-range interacting spin glasses since 
then 0. 

The realistic systems found in nature, however, are 
most likely short-range interacting, quantum SO (3) 
Heisenberg spin systems, which for this reason are es- 
pecially appealing from the physical point of view. In- 
terestingly, however, after more than thirty years, very 
little is known about the properties of short-range quan- 
tum spin-glasses, especially with SO(3) symmetry. Nu- 
merical calculations exist Q , but very few analytical ap- 
proaches are available. Interesting results, however, have 
been obtained in related systems: a quantum, long-range 
SG model [1] has been solved recently @, whereas a 
Landau-Ginzburg, phenomenological approach has been 
developed for quantum, short-range interacting, rotor 
and Ising SG models [10| . 

Nevertheless, still a number of open questions concern- 



ing quantum short-range SO (3) SG systems remain: is a 
mean-field approach possible; is there a SG phase; is it 
stable; is it replica symmetric; is there any basic clash 
between replica symmetry and stability? 

It is natural to expect that a continuum description 
would be useful for a short-range SG system, since the 
nearest neighbor interactions become just derivatives in 
this limit. A serious obstacle arises, however when taking 
the continuum limit of a quantum system. This is con- 
nected with the quantum Berry phases, which in general 
would not cancel when summed over the lattice [TO, [H| • 

In this work we consider a model for a disordered, 
short-range SO (3) quantum spin system in which we cir- 
cumvent this problem and obtain a continuum model 
that will enable us to address the previous questions. 
This is achieved by introducing disorder as a perturbation 
of an antiferromagnetic (AF) 2D Heisenberg model, for 
which the sum of the quantum phases is known to cancel 
[l2l [l3j . The situation is completely different from the 
original EA model, where J = 0, and consequently the 
disorder cannot be taken as a perturbation of a Heisen- 
berg system Q. Using the continuum description, we 
extract the T x J phase diagram of the system from 
the mean field solution, which presents replica symmetry. 
This exhibits a Neel phase at T = 0, whose quantum crit- 
ical point is displaced from its original value in the pure 
system. It also contains spin glass (SG) and paramag- 
netic (PM) phases, whose thermodynamical stability is 
demonstrated by a careful analysis of the Hessian matrix 
of the average free-energy. 

There is in addition an appealing physical motivation 
for this model, in connection to the high-Tc cuprates. 
Indeed, these materials, when undoped, are 2D Heisen- 
berg antiferromagnets, which upon doping, develop a SG 
phase before becoming superconductors. Our model, de- 
scribing precisely the AF-SG transition, is therefore po- 
tentially useful for studying the magnetic fluctuations of 
such materials. 

The model consists of an SO (3) quantum Heisenberg- 
like hamiltonian, containing only nearest neighbor inter- 
actions of the spin operators Si, on the sites of a 2D 
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square lattice of spacing a. The couplings Jy are random 
and associated with a Gaussian probability distribution 
P[Jij] with variance AJ and centered in J > 0, such 
that AJ <C J. We consider the quenched situation, in 
which, according to the replica method [3j the average 
free-energy is given by F = -k B Tlim n — >0 ~ an ~ 1), 
where Z n is the replicated partition function for a given 
configuration of couplings Jij and [Z n ] av is the average 
thereof with the Gaussian distribution. 

Using the coherent spin states |fif(r)), such that 
(Of (r)|Sf \ftf (t)} = STJf (r) (i: lattice sites, a: repli- 
cas, r: euclidian time, S: spin quantum number) [l3j |. we 
may express Z n as a functional integral over the classical 
spin n"(r). The average over the disordered couplings 
Jij can then be performed, yielding 



[Z n ] 



.in 



L j a d,T 



(1) 



where Lj >A = £ ijQ ^ „ + S' 2 J£ (r) ■ Slf (r) + L A . 

In this expression, Lf a are the quantum phases [14j | 
and La is a quartic interaction term, proportional to 
(AJ) 2 , which is generated by the Gaussian average. We 
use the standard Hubbard-Stratonovitch procedure in or- 
der to replace the quartic interaction by a trilinear inter- 
action of flf a (r) with the variables Qf^ b {T, r'), where ab 
are SO (3) indices. 

This is no longer a disordered system. The disorder, 
which was originally present manifests now through the 
interaction term, proportional to (AJ) 2 . In the absence 
of disorder, we would have A J — > and [Z n ] av would 
reduce to the usual coherent spin representation of the 
AF Heisenberg model, with a coupling J > 0, [13, EH • 

Since we are only considering the weakly disordered 
case (AJ <c J) our model, described by the effective 
lagrangian in (TTJ) is a perturbation of the AF 2D quan- 
tum Heisenberg model. This means we can decompose 
the classical spin f2" a (r) into antiferromagnetic and fer- 
romagnetic fluctuations as in that model [14j . Using this, 
then it follows that the sum of the quantum Berry phases, 
Lf a , over all the lattice sites vanishes, as in the pure sys- 
tem [Hill. 

We can therefore take the continuum limit in the usual 
way as in the pure AF 2D quantum Heisenberg model 
HU EH EH obtaining an SO (3) generalized relativistic 
nonlinear sigma model (NLSM). This contains the field 
?i Q = (<7 Q , 7? Q ), which is the continuum limit of the (stag- 
gered) spin £l a and satisfies the constraint n a ■ ri a = p s ; 
where p s = S 2 J. The generalized NLSM also con- 
tains a trilinear interaction of ri a with the Hubbard- 
Stratonovitch field (t, t'), which is proportional to 
(A J) 2 and corresponds to L A . 

Notice that a null value for J, as we have in the EA 
model [|[ would make the perturbation around a NLSM 
meaningless. A negative value, on the other hand, would 
correspond to the ferromagnetic Heisenberg model, which 
after taking the continuum limit, is associated to the non- 
relativistic NLSM. Here perturbation would be possible, 



however, the Berry's phases would no longer cancel. We 
emphasize, therefore, the enormous difference that ex- 
ists, both from the physical and mathematical points of 
view, in considering J as positive, negative or null in the 
Gaussian distribution of the EA model. 

Using the decomposition (g Q/3 = for a = 0) 
Q a J(r;T, T ') = 8 ab [5^xif;r,r') + q<*?(r; r, r')] and en- 
forcing the NLSM constraint with the lagrangian multi- 
plier field A, as usual, we integrate on the ir a component 
of n a , obtaining 



[Z n U = / VaVxDqVXe 



S e fi [a,x,q,X\ 



(2) 



We evaluate [Z n ] av by means of a stationary phase 
approximation, obtaining 



[Z r ' 



_ „-nS cf f[CT^ x ,m 2 ,9° J f(T-T'),Xox('r-r')] 



(3) 



Cf(r-r') 



where nS c g is the effective action in ([2]) , evaluated at the 
extremant configurations cr"(r, r) = a\ 
X(r,r,r') = Xcx(r - r') and g Q/3 (r,T,r') = 
(we henceforth neglect the "ex" subscript). 

The parameter m 2 , as usual, is a spin gap scale, such 
that its inverse is the correlation length. It will be de- 
termined by the temperature, J and AJ, as we show 
below. The staggered magnetization, a characterizing 
an ordered AF state, is defined as a 1 = lim n _>o - ^2 a <r 2 . 

By taking the limit n — > in ([3]) we immediately realize 
that the average free-energy is given by F = 4<5 e ff- This 
is most conveniently expressed in the space of Matsubara 
frequencies oj r = 27rrT, r e Z. Fourier transforming x's 
and g's we obtain the average free-energy density as a 
functional / = / [a a , to 2 , q a ^(uj r ), x(w r )] , where x( u r) 
and q a P(u) r ) are, respectively, the Fourier components of 
x(t — t') and q a ^{r — r'), for which we use the simplified 
notation \r 



x(cj r ) and q"@ = q a/3 (u) r ). 



From ([T]), we can show that Q" p (r, r') = 
(Sf (t)S![ (t')) and therefore, according to the previous 
decomposition of Q into x's and q's, we can identify xo 
as the static magnetic susceptibility, whereas the inte- 
grated susceptibility is given by \i = Xr- The EA 
order parameter for the SG phase [l], l3j , accordingly, is 
given by q EA = Tq , where q = lim„^ „ (w 1 _ 1) E Q/ 3 % P ■ 

By taking the variations of / with respect to the vari- 
ables a a , to 2 , q" 13 , Xr, we obtain the mean field equations 
(MFE). These possess a replica symmetric solution that 
will allow us to determine the phase diagram of the sys- 
tem. 

Let us begin with the search for an ordered Neel phase. 
The MFE imply that this may only occur at T = 0, in 
agreement with (l6j . Indeed, we find that on the line 
(T = }Ps > p ): er 2 = \\p s - p ], qf| = ^-[p s - p Q ], 

X £ F ,xt F ^ oo. Here p = ^[l + I[l+i ln(l + 7 )]~ 

where 7 = 3n ysj*) an d A = 1/a. We also find m 2 = 
so the correlation length diverges on this line. 
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The previous results characterize an AF ordered Neel 
phase (a ^ 0, qea 0) 011 t nc unc (T = 0,p s > p ). 
The parameter 7 appears naturally in the calculation. — 
is a measure of the amount of frustration in the system 
and the actual perturbation parameter. Since we are 
working in the regime of weak disorder, we take 7 3> 
1. We see again that a disorder perturbation would be 
impossible in the original EA model, where 7 = 0. In 
the absence of disorder (A J = 0, 7 — > 00), p — * po(0) = 
which is the well-known quantum critical coupling 
determining the boundary of the AF phase in the pure 
2D AF Heisenberg model at T = 0, 111 El . The effect 
of disorder on the AF phase is to displace the quantum 
critical point (QCP) to the right. This result should be 
expected on physical grounds: in the presence of disorder 
a larger coupling is required, to stabilize an ordered AF 
phase. 

We now search for PM and SG phases. In both of 
them we have a = 0. From the MFE we may determine 
Xr and q r . We find, in particular, go — 0, for to 2 > 
and qo — (3/A)[hIq 



to 2 , 



2 A 

expressions, to,q = — 
Since 5ea = Tq~o it follows that 



m 2 ] > 0, for to 2 < to 2 ,. In these 
ln(l+7)] and A = 7 / (6irp s A 2 ). 

q~o is also a SG order 
parameter and we conclude that the former phase (m 2 > 
to 2 ,) is paramagnetic (a — 0, qEA = 0), whereas the latter 
(to 2 < to 2 ,) is a SG phase (a = 0, (?ea 7^ 0). The phase 
transition occurs at to 2 = m§. In the unperturbed limit 
where the disorder is removed (7 — > 00), we would have 
to,q = and the SG phase would no longer exist. 

We can determine the susceptibilities from the MFE. 



and 



> to 2 ,), we get xf M = 3^ 
3^ - Y(to 2 ,T,p s ,7) = A(to 2 ,T,p s ,7), where 
the function Y(to 2 , T, p s , 7) is explicitly obtained. It has 
the properties Y T ^> A Q an d y J_ [p / p s ] an( j a l so 



In the PM phase (to 
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ln(l + 7), implying that the critical 



X P o M 



satisfies the Curie law 



Ml 



at high- 

- Po/Ps], 



temperatures and diverges as Xo ~ ' 3T l 
for T — > 0. This is the expected behavior for (p s > p ), 
where the AF phase appears. For (p s < p ), conversely, 
we will see that the PM-SG phase transition occurs at a 
finite T c (Fig.l) and the previous expression is no longer 
valid. 

We now turn to the the SG phase (to 2 < to,q). Using 
the MFE we obtain Xo G = A - AAm and xf G = w ~ 
AAmo, where Atoo = [to-q — to 2 ] . 

We can identify a clear cusp at the transition both 
in the integrated and static susceptibilities. This is an 
important result, since the presence of these cusps is a 
benchmark of the SG transition and has been experimen- 
tally observed in many materials presenting a SG phase 

a. 

We may determine the critical curve T c x p s by ob- 
serving that the critical condition to 2 = m 2 , implies 
Xcr - X{ml,T c ,p s ,i) = 0. For T c < A, which corre- 
sponds the situation found in realistic systems, this be- 



comes, near the quantum critical point (p s < po) 



Tc 
2tt 



A 



ln(-) -ln(l + 7 ) 



= p Q -p s . (4) 
We plot the corresponding phase diagram in Fig.l. 
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FIG. 1: Phase diagram for 7 = 10 2 ,A = 10 3 . The critical 
curve corresponds to @ and is valid near the QCP po (solid 
curve). po(0) = A/2ir is the QCP of the pure AF system. 
Notice that disorder besides creating the SG phase, displaces 
the QCP to the right. The value ascribed to A is a realistic 
one in K (A — > ^-A; v s : spin- wave velocity). The result- 
ing temperatures naturally appear with the correct order of 
magnitude, in K, found in real SG systems [J. 

The critical behavior of relevant quantities may be de- 
termined by analyzing the function Y(m,Q,T, p s ,j) for 



T ~ T c and to 2 
for p s < po, 



This yields, near the transition, 



X 



Yc 



47rA 



T - T c 
T r . 



[Po 



From this we can fully determine the critical behavior of 
the SG order parameter and susceptibilities (see Fig. 2). 
For p s > po and T > 0, we always have to 2 > to 2 , and 
a — <Zo — i.e. the system is in the PM phase for any 
finite temperature. A quantum phase transition occurs 
at the point (T — 0,p s — po), connecting the SG phase 
to the AF phase. 

Let us now examine the AJ dependence of the phase 
diagram. For this we make AJ — > AJ(1 + e), with |e| -c 
1, for fixed p s , and study how the relevant quantities 
change. We find, for p s < po, 



Ap e - Ap Tl - T c 



A? 



Atoq 



A 



Ap 



T c 



■In 7, 



47rApo 27T7Ap 

(6) 

where Ap = p - p s . Also, [^ G ] e - ^ G = ^- In 7. We 
see that increasing the amount of disorder, through an 
increment of the Gaussian width, will increase po and, for 
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FIG. 2: The static susceptibility for different values of p s : 
p s /po(0) =(a)1.005 ,(b)0.993, (c)0.981. The cusps, character- 
istic of the SG transition, occur at the corresponding critical 
temperatures (in K). xo is in 

a fixed p g , also T c , m 2 , and q^ G . Conversely, decreasing 
the amount of disorder by narrowing the Gaussian width 
will produce the opposite effects. 

We finally consider the question of the thermodynamic 
stability of the phases. We will focus on the SG and PM 
phases. The stability of the AF phase should not be a 
problem and will be considered elsewhere, in an extended 
version of this paper. The stability of the SG phase is the 
main concern here. For investigating this point, we have 
determined the Hessian matrix of the free-energy density 
/K; fP(wo), ...g^K)...; m 2 ; x (w ), -,xK), ...]. 

This is a matrix with entries of dimensions 
lim„^ [«; n(n— l)o, ...,n(n — l) r , ...; 1; lo, l r , ■•■] corre- 
sponding, respectively, to derivatives with respect to each 
of the above variables. A sufficient condition for the mean 
field solution to be a local minimum is to have all the 
principal minors of the Hessian positive. This would rule 



out the usual instabilities found in long-range replica- 
symmetric solutions, but of course, not meta-stability. 
This would deserve further investigation. 

We have carefully evaluated each of these determinants 
in the limit n —> 0, for a = (PM and SG phases). 
We obtain D a = 1. All the remaining principal minors, 
namely, D qo ,..., D Qr ,..., D m 2,D Xo ,..., D Xr ,... can be writ- 
ten in the form 

£(7-G )+TO), (7) 

where £ and rj are positive real numbers (in the case 
of D qo , for instance, we have £ = r\ = I/7) and G is 
such that in the PM phase Go < 7 and in the SG phase 
G = 7. 

In the PM phase we have qo = and Go < 7, therefore 
it follows that all the principal minors are positive. In 
the SG phase, conversely, (7 — Go) = and q > and 
we conclude that also in the SG phase all the principal 
minors are positive. At the transition, all of them, except 
D(j , vanish, since both (7 — Go) = and % = 0. The 
above result establishes the thermodynamic stability of 
the SG and PM phases presented above. 

Our results, based on a mean-field approach to a short- 
range interacting, weakly disordered, SO (3) quantum 
spin system, clearly show the existence of a stable SG 
phase, at a finite T. The behavior of the susceptibili- 
ties, exhibiting the characteristic cusps at the transition 
is an evident manifestation of it. The solution has all 
q a @ equal, being therefore replica symmetric. The fact 
that it is stable seems to indicate that, in the case of 
short-ranged interactions, there is no basic clash between 
replica symmetry and the stability of the mean-field so- 
lution. 
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